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Abstract 

We consider the goodness of fit testing problem for stochastic dif- 
ferential equation with small diffusion coefficient. The basic hypoth- 
esis is always simple and it is described by the known trend coef- 
ficient. We propose several tests of the type of Cramer-von Mises, 
Kolmogorov-Smirnov and Chi-Square. The power functions of these 
tests we study for a special classes of close alternatives. We discuss 
the construction of the goodness of fit test based on the local time and 
the possibility of the construction of asymptotically distribution free 
tests in the case of composite basic hypothesis. 
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1 Introduction 

We consider the construction of the goodness-of-fit (GoF) tests for dynamical 
system with small noise, i.e., the observations = {Xt, <t < T} are from 
the homogeneous stochastic differential equation 

dXt = S (Xt) dt + sa (Xt) dWt, Xo = xo, 0<t<T (1) 

with deterministic initial value xq and known diffusion coefficient e'^a (■)^ > 0. 
All statistical inference concerns the trend coefficient S {■) only. We have two 
hypotheses: the basic hypothesis in our consideration is always simple 

^n: S(-) = SJ-) 
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and the alternative corresponds to the process ([T]) with a different trend 
coefficient S {■) ^ Sq (■). As usual in GoF testing there are two problems. 
The first one is to find the threshold which provides the asymptotic size a 
of the test and the second is to describe the behavior of the power function 
for some classes of alternatives. There are different ways to present such 
(nonparametric) alternatives and we discuss below the choices of alternatives. 
The problem considered corresponds to the asymptotics of the small noise, 
i.e., we study the properties of the tests as e — * 0. We suppose that the trend 
and diffusion coefficients satisfy the Lipshits condition 

\So{x)~So{y)\ + W{x)-cr{y)\<L \x - y\ (2) 

hence the equation ([T]) has a unique strong solution [T3] and we have the 
estimates 



sup \Xt-xt\> 6} <Ce^^, sup Eo\Xt-Xt\ <Ce\ (3) 

0<t<T J 0<t<T 

where C > 0, c > are some generic constants (see [5], ^Q\). Here the func- 
tion {xt, <t < T} is solution of the limit (deterministic) ordinary equation 

d T+ 

-^ = So{xt), xo, 0<t<T. (4) 

Our goal is to present some goodness-of-fit tests for this stochastic model 
which are similar to the well-known in classical statistics Cramer-von Mises 
(C-vM), Kolmogorov-Smirnov (K-S) and Chi-Square (Ch-S) tests. Remind 
that these classical tests are distribution free, i.e., their limit distributions do 
not depend on the basic hypothesis and therefore the problem of the choice 
of the threshold is universal for all tests of such types. Moreover these tests 
are consistent against any fixed alternative. The GoF tests proposed below 
for the model ([T]) have the similar properties. 

Let us recall the basic properties of the classical tests. Suppose that we ob- 
serve n independent identically distributed random variables (Xi, . . . , Xn) = 
X" with continuous distribution function F (x) and the basic hypothesis is 
simple : 

J^o- F{x) = Fo{x), X e M. 

Then the Cramer-von Mises and Kolmogorov-Smirnov Dn statistics are 

/°° r - 1 2 

F„ (x) - Fo {x) dFo (x) , D„ = sup F„ (x) - Fq (x) 

respectively. Here 

1 " 

3=1 
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is the empirical distribution function. Let us denote by {Wq (s) , < s < 1} 
a Brownian bridge, i.e., a continuous Gaussian process with 

EWo (s) = 0, EWo (s) Wo{t)=tAs- st. 

Then the hmit behavior of these statistics can be described with the help of 
this process as follows 



Dr. 



Hence the corresponding C-vM and K-S tests 



sup \Wq (s) 

0<s<l 



with constants Cq, da defined by the equations 
»i 



Wo isY ds > c. 



a. 



sup \Wq {s)\ > do 

0<s<l 



a 



are of asymptotic size a. We see that these tests are distribution-free (the 
limit distributions do not depend on the function F^.. (•)) and are consistent 
against any fixed alternative (see, for example, Durbin [3], Lehmann and 
Romano [Hj). 

It is interesting to study these tests for non degenerate set of alternatives, 
i.e., for alternatives with limit power function is grater than a and less than 
1. It can be realized on the close nonparametric alternatives of the special 
form making this problem asymptotically equivalent to the signal in Gaussian 
noise problem. Let us put 



F{x) = Fo{x) + ^ [ 



h{Fo{y)) dFoiy), 



where the function h{-) describes the alternatives. We suppose that 



his) ds = 0, 



h (s) ds < oo. 



Then we have the following convergence (under fixed alternative, given by 
the function /i (■)): 



Uo 



sup 

0<s<l 



h (v) dv + Wo (s) 

s 

h (f ) df + Wo (s 



ds, 
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We see that this problem is asymptotically equivalent to the following signal 
in Gaussian noise problem: 



dYs = K (s) ds + dWo (s) , < s < 1. 



(5) 



Indeed, if we use the statistics 



D = sup \Ys 



Jo 



0<s<l 



then under hypothesis h{-) = and alternative h{-) ^ the distributions 
of these statistics coincide with the limit distributions of and Z)„ under 
hypothesis and alternative respectively. 



2.1 Choice of the thresholds 

We test the basic simple hypothesis J^o and our goal is to study the GoF 
tests of asymptotic size a G (0, 1). Let us denote the class of such tests as 



where Eq is expectation under hypothesis. We use the following regularity 
condition. 

Condition TZ. The function Sq (■) has two continuous bounded deriva- 
tives S'q (x) , Sq (x) and the function a (x) has one continuous bounded deriva- 
tive a' (x) and the both functions are positive : Sq (x) > 0, (x) > for 
X > Xo . 

Remind that under this condition the equation ([T]) has a unique strong 
solution and this solution converges uniformly on t G [0, T] to the determinis- 
tic solution of the equation (jll) (see, e.g., [15] [IDj)- Moreover, the stochastic 
process Xt is differentiable w.r.t. e at the point e = and its derivative x|^^ 
satisfies the linear equation 



For the proof see, e.g., [lU], Lemma 3.3. 

To construct the C-vM and K-S type tests we use the statistics 



2 C-vM and K-S type tests 



/C„ = {7/>, : Eoij,{X') = a + o{l)} 



dx[^^ = S'o (xt) xj^^ dt + a (xt) dWt, x[,^^ = 0. 



(6) 




2 



4 



and 



Is = 



.5'o {xt) 



-1/2 






sup 

0<t<T 


eSo{xt) 



respectively. 

Below Ca and ba are solutions of the equat 



ions 



sup \wy\ > he 

0<D<1 



P |y wl(lv> Coj = a, P 
where w^,, < f < 1 is some Wiener process. 

Proposition 1 Let the condition TZ he fulfilled then the tests i/j^ = 1 
and (h^ = 1 belong to the class /Cq. 



L , , , belong to the class K, 

Proof. The stochastiLv piw^^^oo c y^'^t — ^t) 
formly on t G [0, T] to the limit x^^^ - solution of 
for any k > 



(7) 



process e ^ {Xt — Xt) converges in probability 
,„ *the linear equation ([6]), i.e., 



uni- 



Xt-xt _ (1) 
■^t 



Po { sup 

[0<t<T 

This solution can be written explicitly 



> K 



0. 



as 



x\ = 



exp 



(x^) dv [ cr (xs) dWs. 
Jo KJs ) 

Below we follow [TU], where the similar calculus were done. Using ([21) we can 
write 

»S'q (x^,) dx^ 



S'q {xy) dv 



So (xy) dv 



dv 



o.x^ — 

5 '-'0 \^v) 



: r^^dx= r (inSoix))' dx = \n 



So (Xs, 



Hence 



X 



(1) 



Sq {Xs 



ds 



where W {■) is some Wiener process. Further 

i-T 



cr [Xsj dt 



\eSQ{xt) 



W 



r^ixsf 

lo So{Xsf 



ds d 



lo So{Xs 



■ds 



: / W {Uf du=ul 

Jo Jo 



wl dv, 
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where we put 



u 



as = V 



So i^s 







5*0 {xs 



ds = V ut- 



Here = Uj, ' W {utv) , < f < 1 is a Wiener process. 
Hence (under hypothesis Mq) 



dv. 



For the statistic 7^ the similar consideration leads to the relation 



7£ — > sup \w^\ 

0<v<l 



Therefore for the first type errors we have 



and 



ae (S,) = Po {5, > Ca} — ^ P <! ^ wldv> 



ae (7e) = Po {le > ba} > P <J SUp \Wy\ > bo 

n<v<i 



a, 



a. 



The Proposition [T] is proved. 
2.2 Similar tests. 

Note that the similar result we have if we use the statistics 



\So{Xt 



dt 



\eSo{X, 



and 



7e 



\So {Xt 



dt 



-1/2 




Xt - xt 


sup 

0<t<T 


eSoiXt) 



because, as we mentioned above, the process Xt converges uniformly on t G 
[0, T] to the deterministic solution Xt and this implies the convergence 



5. 



wl dv, 



le 



sup \Wv\ 
0<v<l 
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These tests can be slightly simplified, if we replace the equation (jlj) by 
the following one 

dX 

' SoiXt), Xo = Xo, (9) 



dt 



Xt = xo+ So (Xs) ds, 0<t<T, 



i.e., we put 



(we need not to solve the equation (jlj), just to calculate the integral, using 
observations) and introduce the statistic 

2 



1 -2 



cr{Xtfdt 



\2 / Xt — Xt 



dt 



Then under hypothesis 



e-^ / a{Xtf(Xt-Xt) dt= / a(X, 



(y {Xs) dW, 



dt 



W[ / a{Xsfds 



-I 2 



d / a{X,fds 



TT 



W (rY dr 



W {rf dr = {T°f 



wl df , 



where 



Hence we have the convergence 



tt = o [Xsf ds — > = / cr {xsf' ds 



w„ dv. 



Of course, we have the distribution free limit for the corresponding statistic 
% too. 



2.3 Partially observed linear system 

Suppose that we observe a random process X'^ = {Xt, < t < T} and we 
have to test the hypothesis J^q that this process comes from the following 
linear partially observed system 

dYt = AYtdt + eBtdVt, Yo = yo, 0<t<T, 
dXt = CtYtdt + eatdWt, Xq = 0, < t < T, 
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where At,Bt, Ct and at are known functions and Vt and Wt are independent 
Wiener processes. We suppose as well that these functions satisfy the usual 
conditions which allows us to write the equations of filtration (see Liptser 
and Shiryaev [15]). 

If the hypothesis is true then according to well known Kalman-Bucy the- 
ory (see Liptser and Shiryaev pjj, Theorem 10.1) the conditional expectation 

= Eq (1^1X3, < s < t) satisfies the equations 

dMi = AtMtdt + [dXt - CtMtdt] , 



dt 



with initial values Mq = EqFo = Vo and 7 (0) = Eq (Iq - EqIq)^ = 0. Note 
that if we put = 7 (t) e~^, then this system can be rewritten as 



dM. = At Mt dt 



CtTt 
i2 r2 



|dX, -C,A/,di], M„ = y„, 



dt 



at 



t 1 



0. 



Remind as well that the observed process admits the representation 

dXt = Ct Mt dt + eat dWt, 0<t<T, 

where Wt is innovation Wiener process defined by this equality. This repre- 
sentation suggests to define the statistic 
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e I a{t) dt 




T 



a{ty 



Xt 



M, ds 



dt. 



Elementary calculations yield 



a{tf 



Xt - / C.MAs 



a{ty 



w 



a. 



'ds 



dt 
2 



a{tf 



a, dW, 



dt 



dt= (^J^ a {tf dt^ wldv. 



Hence the statistic 5^ (under hypothesis) is 

•1 



5. 



w„ dv 



and the corresponding test ipe 



is distribution free. 



2.4 Alternatives 

Let us consider nonparametric alternatives S {■) ^ 5'o(-), which correspond 
to the equation 

dXt = S {Xt) dt + ea {Xt) Xo = xq, 0<t<T. 

We suppose that the functions S {■) and cr (■) satisfy the conditions, 

IS* (x) - S{y)\ + \a (x) - cr {y)\ < L\x-y\, 
\S{x)\ + \a{x)\ < L(l + 

The set of such functions we denote as S (L) and the hmit of the stochastic 
process Xf, < t < T we write as Xt (S) ,0 < t < T. 

There are many ways to introduce the class of alternatives. We consider 
two of them. Let us define the sets 

= {S {■) e S (L) : ||a;. (5) -x.ll > r}, 
gr = {Si-)eSiL): \\Six,)-Soix,)\\>r}, 

where ||-|| is £2 (0, T)-norm, say, 

\\S (x,) - So {x,)f = f [S {xt) - So {xt)f dt. 
Jo 

Here Xt, < < t < T is solution of the equation (jlj) (under hypothesis J^). 
Let us start with the problem 

Jfo '■ S {■) = So (■) , 

J^i : S {■) E Tr 

To show the consistency of the test ■j/'g note that > ^^s^^ ||X. — with 
some K > 0. We have 

inf Vs {be > c„} > inf Ps {^e-^ \\X, -x,\\> 

> inf Ps {e-^ \\x, (S) - x,\\ - e'^ \\X, - x, (S) \\ > ^} 

> 1 - sup Ps {e-^ \\X, - X, (S) II > e-^ \\x, (S) -x,\\- K-^^/^} 

> 1 _ (e- V - n-^J^) sup Es r ( ^^l^^li^) ' dt > 1 - Ce\ 

^ ' SGTr Jo \ £ J 

Therefore even if r = — > such that e'^r^ —>■ 00, then the C-vM 
test is minimax consistent, i.e., its power function tends to 1 uniformly on 
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S {■) e J-'r- The same is true for the Kolmogorov-Smirnov test. It is not so 
in the case of the hypotheses testing problem 

J^o : ^(■) = ^o(-), 

J^i : S {■) E Qr 

It will be more convenient to write S (x) — Sq (x) as eh (x) a [x)^ Sq {x)~^, 
where /i (•) is such that S {•) G Qr- This corresponds to the model 

dXt ^ So (Xt) dt + e ^ ^^c^!"}^'^ + ea (Xt) dWt, 0<t<T (10) 

with the same initial value Xq = Xq. The case h{-) = corresponds to the 
hypothesis Mq. We start with the study of the limit behavior of the statistic 
5e under a fixed alternative h{-). 

Let us denote by the solution of the equation 

and write 

Xt — Xt Xt — x'/ .X'/ — Xt 



It is easy to see that 



— ^ — >xf\ 0<t<T 



e 

in probability and the direct calculations yield 

x1 -xt_ /■* So {x':) -Soix,) , a {x':) ' h 



Jo £ Jo 



Hence 2t— ^ converges to the function x'l which is solution of the equation 
Therefore 

x1 = exp 1 1' S', (x,) dv^ h (xs) ds 

= -^0 (xt) I |-7^2 ^ i^s) ds 

Jo bo [Xs) 
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and 



eSo {xt) 



So (xt) 



So {XsY 



h{xs) ds + wl [ ^ 
\Jo So 

Jxo So {xY V Jxo Sq 



[X. 



[X 



ds 



da; 



Let us put 



u {x) 



dy 



[0,ut] 



'xo SoivY 

and denote x (u) the function inverse to u (x) These relations yield the fol- 
lowing representation for the limit of the test statistic 



Jo 



Z{xtfa{xty 
So {xtf 



dt 



r 

Jo 



r — 

,Jxo So 



XT 



XO 



Z(x) 



Un 



ut 



xo 



So {yy 



dy 



— Un 



h{x{z))dz + W{u) 



du 



K (S) ds + Wy 



ds, 



(11) 



1 /2 

where we put u = utv and denoted /i* (s) = Uj, h{x{uTs)) and w. 



u 



-1/2 



T 



W {utv). It is easy to see that w^jO < v < 1 is standard Wiener 
process. 

Remind that if the observed process is of type signal in white Gaussian 
noise: 

dYs^h(s) ds + dws, < s < 1, 

then the natural distance between hypothesis h{s) = and h{s) ^ is 
(0,1): 



/ 

Jo 



h{sf ds 



In our case it corresponds to 



fK{sf ds^ rhixf^. 

Jo Jxo So 



n dX > P^. 

{xf - 
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The convergence ffTT]) provides the hmit of the power function given in 
the next proposition. 

Proposition 2 let the condition TZ he fulfilled, then for any function h (■) G 
Tir we have 



/l* (S) dS + Wy 



ds> cA +o{l). 



Let us define the composite (nonparametric) ahernative as 
M[: h{-)enr, 

where 



hi-] 



h (x.) a {x,y 



So (x.) 



> ^ 



Note that x, here is "from hypothesis J^". As h (■) is an arbitrary function, 
this alternative coincides with Qr- Let us show that without regularity con- 
ditions the problem of hypotheses testing is degenerate in the following sens: 
we have for the power of the test 



inf P (4, h) 



a 



as e — s> even for fixed r > 0. The condition S {■) G S (L) already provides 
the control of the first derivative, that is why we have to weaken it slightly. 
At particularly, we suppose that S {■) G 5 (L), but L > can be not the 
same for different functions S {■). 
Let us put 



So (x) 



aix 



cos {n (x — Xq)) 



then, 



\Sn (x) — S'o (x) I = ec \So (x) cos {n (x — xo))| 

and L = Keen with some K > 0. Further, using 2 cos^ if = 1 
obtain 



cos {2(f) we 



hn [x) — ^ dx 



So {xY 



c 
~2 



XT 



So (a;) 



XT 



— \ COS {fin (x — Xo)) 



a [xY 
So [x) 



dx 



cr X 



2 dx > — 



XT 



So (x) 



dx 



Xo 



a X 
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as n ^ oo. Remind that 5*0 (x) and o (x) are continuous positive functions. 
The constant c = c (e) > can be chosen from the condition 




Hence /?.„ (•) G Tir- On the other hand 



uniformly in x G [xq, xt\ if we put, say, n = n{e) = c{ef . Therefore 








lim inf 3 (5^, K) < lim inf 3 (5^, h„) = a. 



Hence if we use the introduced above statistics then this hypotheses testing 
problem will be asymptotically degenerated. It can be shown that this is not 
the particular property of these two tests, but is true for any other tests too. 



then it is natural to replace Xt (S*) by its estimate Xt and this leads to the test 
= l{||x.-x.||>ec} similar to the introduced above C-vM type test. Remind 
that it was shown before that Xt is the best in different senses estimator of 
Xf (S*) (see [lO], Section 4.3). Therefore we can think that if the alternative 
is defined by the relation 



then for construction of the good test statistics we have to replace S (xt) by 
some estimator. Forget for instant that the Wiener process is not differen- 
tiable and write the observed process as 



3 Chi Square Test 



If the alternative is defined by the inequality 



||x. (S) — x.ll > r. 



S (x.) - 5*0 (x.) 
ct(x.) 



> r. 



Xt = S (Xt) + ea (Xt) Wt 
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Then (formally!) Xt ^ S {xf) as e — > and we can consider Xf as "esti- 
mator" of S (xt). Let {(fj (•) , J = 0, ±1, ±2 . . .} form an orthonormal base in 
£2 (0, T), then by Parseval identity (formally!) we have the following equality 
for the corresponding (modified 5*0 (xt) Sq (Xt)) statistic 



— So {Xt 



where the Fourier coefficients yj^e we can write as follows 



£(^{Xt] 



Xt-So{Xt) 



I 

Jo 



ea {Xt 



[dXt-So{Xt)dt] 



This last integral has mathematical meaning and starting from this definition 
of Uj^s we can introduce the statistic 



\j\<m 

Note that under the random variables yj,£,j = 0, ±1, ±2 . . . are indepen- 
dent Gaussian 



yj,- 



ifj {t) dWt 



M{0,1). 



Therefore the statistic [Vj^^ — l] has Chi-Square distribution and the 

equation for 

Po {6* > Ca} = a 

can be easily solved. Moreover, m — > 00 we have the convergence (under 
hypothesis) 



\j\<m 



Hence we can introduce the Ch-S test as 

where m = m (e) — >• 00 and Za is 1 — a quantile of the Gaussian J\f (0, 1) law. 
This leads us to the following result. 

Proposition 3 Let us suppose that a (x)^ > 0,x & R, then the test ijj* {X^) 
belongs to K^- 
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We see that we need not even to use the convergence Xf to Xt and the 
proposition is vahd, say, for ergodic diffusion processes with e — 1 and the 
asymptotic T ^ oo. The choice of asymptotic is important in the calculation 
of the power function. 

Suppose that the observed process has the trend coefficient S (•) ^ So (•) 
and the condition TZ is fulfilled, then 



£Cr[^t) Jo 



where 



E4= 



T 



S{Xt)-So{Xt 



dt 



S (x ) — Sn (x 



a{x.) 



Hence for any fixed contiguous alternative S {x) — S (x) + eh (x) a {x) we 
have 



\3\<ra 



|j|<m 



|j|<m 



As m — > oo the relation 



|j|<m 



\\h{x)f{\ + o{\))>-{\ + o{\)) 



holds. To have non degenerate limit 

1 ^2 

/4m ^ 



r u 



£2 vim 



:i + o(l)) -^u>l 



\j\<m 

we can put m = r''/ (4£:''). Note that in this case 



E 



'4m 

\j\<m 



2 

r u 



\j\<m 



0. 



Hence for the power function we obtain the limit 

p (s;, h)^p{c>z^-u}, c (0, 1) . 

Another way is to fix first m — oo such that me^^'^ — >■ oo and then to consider 
the alternatives running to the hypothesis : r = r (s) = e\/Am — > 0. 
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Of course, this test is not uniformly consistent with the same exphcation 
as above. 

To have uniformly consistent and minimax GoF testing we need to control 
the derivatives. If we suppose that the function h{-) defining the alternative 
is k times different iable and the £2 norm of the fc-derivative is bounded by 
some constant, then we can show that the test 



|j|<m 



with weights 



Wi = z 1 



m 



2k^ 



2k 



^1/4 



and special choice of m — 00 is asymptotically minimax. The proof is based 
on the approach developed by Ermakov |1] and Ingster and Suslina [9] for 
the model signal in white Gaussian noise. See as well the similar problem for 
Poisson processes studied by Ingster and Kutoyants [8j. 



4 On local time and GoF testing 

The local time of the diffusion process (P) is defined as 

p2 rT 

Arix) =lim— / 1,^ a{Xtf dt 



and admits the Tanaka-Meyer representation (see |18j ) 

r 

At (x) = \Xt — x| — |xo — a;| — / sgn {Xt — x) dXt. 

Jo 

The local time At [x) recently started to play an important role in statistical 
inference [H] , [1] , [12] ■ Note that in ergodic case {e = 1 and T 00), the 
local time is asymptotically normal : 



At (x) 



^Ta (x) 

where / {x) is invariant density and 
djixf = Af{xf E 



fix)] ^M{0,df{xf) 
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Here ^ is random variable with density / {x) and F (x) is its distribution 
function (see [T2], Proposition 1.25). Moreover this normed difference con- 
vergence weakly to the limit Gaussian process in the space of continuous on 
TZ functions vanishing in infinity [I2j, Theorem 4.13. This property can be 
used in the construction of the GoF tests as follows. Let us introduce the 
C-vM and K-S type statistics 

/oo 
riT (x)"^ dx, 7t = sup \r]T (x) \ 
-oo ^ 

and the corresponding tests ipT = 1{5t>cq} 0t = ^{-yrxic}- Then using 
this weak convergence we can define the constants Cq,, da (see [I2], section 5.4 
and Gassem [6]). 

We can consider the similar problem in the asymptotics of small noise, 
i.e., to use the limit behavior of the local time in the construction of the GoF 
tests. 

Let us introduce the space {xq, xt) of square integrable functions on 
[xc^t], where Xt is solution of the ordinary differential equation 

dxt _ 

— - — Do \Xt) , Xq, 

at 

where 5*0 (x) > 0. According to (|T2l) we have 



At {x) _ 1 



T 



lim = lim lim / I {\x,-x\<u}(^ {Xtf dt 

= lim— / l{\cc,-x\<u}C^ [Xt) dt=lim— / -— dy 

v^o 2v Jq '- v^o 2u J^^ So [y) 

'^^o2z/J^_^ So{y) So{x) 

t 

We say that the random process rj^ (x) ,xo<x<xt converges weakly 
in {xo,Xt) to the random process r] (x) , Xo < x < Xt if for any function 
h (■) G (xo, Xt) we have 

XX rxT 

h (x) r^g (x) dx =^ / h{x) 1] (x) dx. 

'XO Jxo 

Proposition 4 Let the condition TZ be fulfilled, then for the local time (■) 
we have the weak convergence in C? (xo,xr) 

\2 



^Jxo\S^VV) e^oiy) J \Jxo SoiyY 

'(13) 

where x G [xq, xt] and W {■) is a Wiener process. 
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Proof. Note that the local time allows us to write the equality (see [TH] ) 

r h (Xt) dt = r h (x) 

Jo J-oo e^a(x] 



Suppose that h (■) G C^. We know that 

^0 JO JO So{xt) Sq{x 

and 



io ^ JO 



^0 \Jxo So{y) J 

r h' (x) wlr -^dy ] dx 



Hence if we denote 



e^cr I a; I 



and 



5o(a;) ; 



then we can write 



/i (x) 7]e (x) dx — > I h' (x) ly {g (x)) dx. 



The same time integrating by parts we have 



h (x) r/g (x) dx = / /i' (x) ■j/'^ (x) dx, 



where we put 

ijje (x) = / % iy) dy. 

J X 

Hence 

r-XT 

h' (x) ijji; (x) dx — > / h' (x) W {g (x)) dx. 



Xo 
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We see that the values of h{-) outside of the interval [xq, xt] have no contri- 
bution in the limit. Therefore we have as well the convergence 

XT r^T 

h' {x) il)e{x) dx — > I h'{x)W{g{x)) dx. 

1X0 J XQ 

Remind that this is true for any function h (■) G . Therefore the proposition 
is proved. 

The convergence ( |T3i) suggests the construction of the following test. Let 
Then it can be shown that 

"1 

w'^. df 



and the corresponding test ips = l{5j->c„} is asymptotically distribution free. 

Remark. We see that despite the ergodic case, the local time random 
function has no limit (as process) and for small values of e its behavior is 
close to white noise process. 

Remark. We supposed above that So (x) > for all x G [x^xt]. In the 
case S'o (xq) = the deterministic solution xt = xq and we have the following 
basic hypothesis 

^0 : dXt = ea {Xt) dWt, Xq, < t < T. 

The test can be based on the statistic 

and it is easy to see that (under hypothesis ) 

4 — ^ / wl ds. 



If for some x^. > xq we have S'o (x) > 0, x G [xq, x*) and S'o (x*) = 0, then by 
Lipschitz condition 



dx 1 dx 1 , X* - Xq 



and we see that the equality Xf^ — OC ^ IS impossible (well known property). 



19 



5 On composite basic hypothesis 

Suppose that under hypothesis ^ the observed diffusion process is solution 
of the stochastic differential equation 

dXt = S Xt) dt + ea {Xt) dWt, Xq = xo, < t < T, 

where the trend coefficient 5* (t?, x) is a known function which depends on 
unknown parameter G = 7). Then the limit solution Xt depends on 
the true value i.e., Xt = Xt{'d) and the natural modification of the test 
statistic can be based on the normalized difference 

where as we can take, say, the maximum likelihood estimator. We suppose 
that the functions S (t?, x) and a {x) are positive and sufficiently smooth 
to calculate the derivatives below and to provide the "usual properties of 
estimators". Remind that (under regularity conditions) this estimator is 
consistent and asymptotically normal (see [10], Theorem 2.2). Moreover the 
MLE admits the representation (see [10], Theorem 3.1) 

k_ = % + sir (tf„)- r ,av, + o (.) . 



/o (y{xt{^Q)) 

where do is the true value of d and \t {'do) is the Fisher information: 

It can be shown (see Rabhi [IT]), that in regular (smooth) case the limit 
distribution of ^ 

5e= ! Yt{d,)^dt 

Jo 

coincides with the distribution of the following integral 

" (^0) - It (^o)"^ Xt (^0) T ^^^^Pdiy.^ ' dt. 



a{xs{i^o)) 

where Xt (^9) = -§0 Xt {^) and x^^'^ {d) is solution of the linear equation 

dxl^^ = S' {'&, Xt) xf^ dt + a {xt) <\Wt, x^Q^ = 0, < t < T, 
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The test is no more distribution free, but in some cases it can be done 
asymptotically distribution free (ADF) if the second limit T — cxd is taken 

Another possibility to have an ADF test is to use the estimator process 
"^1,6,0 ^ t < T, where 'dt,e is an estimator constructed by the observations 
X* = {Xg, < s < t}. Of course, we have to suppose that this estimator is 
consistent and asymptotically normal for all values of t G (0, T]. For example, 
in the linear case 

dxt = 'dh (Xt) dt + ea (Xt) dWt, Xq = xo, 0<t<T 

we can take the MLE process 

Sometimes an estimator process can have recurrent structure (see Levanoy 
et al. [TB]). If dt^e is the MLE, then it can be shown (see the similar calculus 
above (El)) the limit of is (below Xs = Xs {'^o)) 

Jo S (i'o, Xs) Jo a [Xs) J 

= / S{^o,xtfZ{t,^ofdt. 
Jo 

The process Z {t, ^q) has stochastic differential 

dZ (t, i9o) = A (t, ^o) dt + B (t, ^o) dWt, Z (0, ^o) = 0. 

with the corresponding random function A (t, t^o) and deterministic B {t, i^o)- 
Hence ^ 

B {t, do)-' (^Z (t, ^o)- j^A {s, do) ds^ = Wt 

and 

* = X ( rs(tf„.x,)B(M„) j * = X 
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Using "empirical versions" of these functions it is possible to construct 
an ADF test based on the following statistics 



V TS{'d,,Xt)B,{t,'&,) 



dt 



w,, dv. 



What is empirical version of stochastic integral we explain below. The inte- 
gral is started at t — e'^, where ji e (0, 1) because for the values t e [0, e] the 
estimator dt^e is not asymptotically normal. 

One else ADF test can be constructed by "compensating" the additional 
random part by the following way. First we rewrite the stochastic integral 
(Ito formula) 



Jo o-{Xs) Jo ^ 



S{^,Xs 



(Xs) 



-dW.. 



'""^ S y) 



dy-e^ 







Jxo o- (y) 

_ S{d,Xs) S{d,Xs) 
Jo 



T CI 



S' X,) a (Xs) - 2S (i?, X,) a' {X,) 



2a {X,) 



dt 



(jiXsY 



ds. 



The last expression for does not contain stochastic integral and we can 

put the estimator, i.e., the random variable (^^e^ is well defined (empirical 
version) . Then introduce the stochastic process 

Note that it can be easily shown that 



sup 

0<t<T 



where 



xti^e) - xt (i?o) ^ 0, Iri^e) ^ It i^o) , H^i'd,) ^ Ho {^o) , 

^ S{'do,Xs (t?o)) 



Ho i^o) 



-dW.. 



/o a {xs i'&o)) 

Hence the stochastic process Y^{t,'&e) converges uniformly on i e [0,7"] to 
the Gaussian process x^^^ (■j?o) and we can use the statistic 

-2 



/ 

Jo 



a{Xt 



S{^e:Xt) 



-dt 



Ut,^sYa{Xty 
s{^e,Xty 



dt 



dv. 
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We see that the test ip^ = 1{4>C(,} based on this statistics is ADF. 

This test has to be consistent against any fixed alternative. Indeed, let 
the observed process be 

dXt = S (Xt) dt + ea (Xt) dWt, Xo = Xq, 0<t<T 

and the limit solution xt {S) satisfies 

g = mf {S)-x, {^)\\ > 0. 

The MLE -t^s in this misspecified situation converges to the value which 
minimizes the Kullback-Leibner distance (see ^Q\, Section 2.6) 

-d^ = arg mf / — — dt. 

Hence 

(t, I) = xi'^ (S) + ^^(^)-^* (^*) + (^^) H, (^,) + o (1) 

and 



■'e 



OO. 



Therefore the test is consistent. 
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